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In this contribution we present the results of a series of investigations of dimensional 
£SJ . reduction, applied to SU(3) gauge theory in 2 + 1 dimensions. We review earlier results, 

, | present a new reduced model with Z3 symmetry, and discuss the results of numerical simil- 

es ' lations of this model. 

§1. Introduction 
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> 



O 



Dimensional reduction is a powerful method to study the long distance behaviour 
of field theories at high temperature. It was first introduced in [1, 2]. It was further 
developed for gauge theories in [3, 4, 5]. The first quantitative studies using lattice 
CN ■ simulations of the reduced model were performed in [6, 7, 8, 9]. Quantitative results 

for QCD were further obtained in [10 - 15]. For a nice review see [16]. 
I In order to find the region of validity of dimensional reduction we have inves- 

tigated in detail the reduction of SU(3) gauge theory in 2+1 dimensions [17, 18, 
19, 20, 21]. In section II we present the method, in Section III some of the results 
are described. In Section IV an alternative form of dimensional reduction is pre- 
sented, which preserves the Z3 symmetry of the action. Section V, finally, gives our 
Q-i, conclusions. 

§2. Perturbative dimensional reduction 

The properties of gauge theories at finite temperature is described by the parti- 
al ' tion function 

' VA a ^Vip^ c Vip^ c e~ s ^'^ (2-1) 

where A^(x) are the gauge fields in the adjoint representation, and where we have 
included fermion fields of flavor / in the fundamental representation of SU(3). The 
action S at vanishing chemical potential is given by 



-l/T 

S = I dx I d ds x 







f 



(2-2) 



The metric is Euclidean and d s is the number of space dimensions, x = (xq,x). The 
physical case is d s = 3. In this contribution we will mainly discuss a simpler theory, 
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where the method of dimensional reduction can be tested in detail, namely pure 
SU(3) gauge theory in 2+1 dimensions (d s = 2) without fermion fields. 

The xo-i n t e g ra ti° n goes from zero to 1/T, where T is the temperature. For 
T — > do, and for distances R 3> 1/T it seems plausible that one can integrate out 
the non static modes of the field perturbatively. It is important to note that the 
perturbation theory of these modes does not have infrared divergences, in contrast 
to the perturbative expansion in the full theory. We thus obtain an effective model 
of the static modes in d s dimensions. Note that the fermion fields, which obey 
Fermi statistics have antiperiodic boundary conditions in xq, and thus have no static 
modes. 

Z = JvA${x)VAf (^e-^WC*)-^)) (2-3) 

It can be shown, that for T large and R ^> 1/T (\p\ -C T), S e jf has a finite number 
of local terms at given order in g 2 /T, \p\/T [5]. The coefficients are determined by 
the renormalized perturbative expansion, and there are thus no new free parameters. 
A pictorial description of dimensional reduction is given in Fig. 1, where the word 
Higgs refers to the field ^4q*. Obviously the reduced theory demands much less 
computertime than the full theory. 
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Fig. 1. Pictorial description of dimensional reduction 
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The static modes are defined by 

A^x) = A$(x) + A™(x); [ 1/T dx A™(x) = 

Jo 



(2-4) 



Note that this splitting is not gauge invariant in general. The physical quantities, 
which we calculate are, however, gauge invariant. Define 



Af(x) = VTMx) 
A s \x) = VT4>(x) 

On the classical (tree) level we set 



S ef f= [ d^x\-[F-(x)} 2 + J2Tr[D. 

J L i=i 



(2-5) 



(2-6) 



This we also call "naive reduction". Because of the ultraviolet behaviour of the 
theory, we should, however, also take quantum effects into account. Restricting to 
the local terms, important for large distance physics, we find to one loop order the 
systematics described in Figure 2. 
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Fig. 2. Contributions to the 2 ,0 4 and 6 terms in the reduced action. In the diagrams, dotted 
lines are static modes, full lines non-static modes. The numerical coefficients can be found in 
the original papers, for 3+1 dimensions in [6, 7, 8, 10, If] and for 2+1 dimensions in [17]. 

Note that in d s = 3, g 2 is dimensionless. For d s = 2, g 2 has dimension energy, and 
4> is dimensionless. In the latter case, obviously higher loop orders are suppressed at 
high T by powers of g 2 /T. 

The form of S e f / calculated through perturbation theory of the non static modes 
in the lattice regularization can be found in [5, 6, 7, 8], and in dimensional regular- 
ization in [10] [11]. Here I present the result for d s = 2, pure gauge theory, which we 
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derived in [16]. We have chosen the static Landau gauge [STALG] 

doAo(x) = 

/ diAi{x)dxQ = 
Jo 

We find to one-loop order in lattice regularization with spacing a 
S 2 ee S eff = J d 2 x + Y J Tr[D^] 2 

L i 

-^f I log 2 - 1 - log aT]Tr4? + |W 



(2-7) 



(2-8) 



The effective model is a Higgs gauge model in two dimensions, with the Higgs field 
(f) in the adjoint representation. The action has a particular symmetry, i?-symmetry, 
under <fi — ► — 0, which comes from the T-reflection symmetry in the 2+1 dimensional 
theory. Including fermions in the original model will only make a change of the value 
of the coefficients. At finite density, however, i?-symmetry and hermiticity is broken, 
and there are also terms with odd powers of 4> and imaginary coefficients [15]. 

There is a logarithmic divergence in the coefficient of the quadratic term in 4>. 
This will be cancelled by the ultraviolet divergence of the two dimensional model. 
The perturbation theory of the model is infrared divergent. And although the pure 
gauge theory in two dimensions has been solved on the lattice [22], and in the contin- 
uum in the large N limit [23] , there is no known solution for 2d gauge theory coupled 
to Higgs fields in the adjoint representation. Instead we solve the model numerically, 
using lattice Monte Carlo simulations. As a lattice version of this action we choose 

S 2L = p 3 LoY / (l-^ReTrU p ] 



(2-9) 



7tL /?3 



-log2- 1 + logLo 



Tr<p 2 (x) + 



8tt/?; 



The parameters (3$ and Lq which appear in this action are related to the coupling 
g\ and the temperature T of QCD2+1 by 



/3 3 = — 5 L = l/aT 
a 9l 



(2-10) 



where a is the lattice spacing. 

The results for physical quantities using this effective action are compared with 
the full 2+1 dimensional SU(3) gauge theory, calculated with the lattice action 



S3L = /? 3 £E' i 

x ii<v 



1 



ReTrU p (x, fj,, v) 



(2-11) 
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with the same spacing. 

The matrices U p and U p are the products oiU(x, i) resp U(x, i) around a plaque- 
tte in 3 resp 2 dimensions, U(x, i) and U(x, i) £ SU(3). The matrix <j>(x) is, however, 
in the algebra of SU(3); thus the ^-symmetry of S%l is broken in S^L- 

We also define a dimensionless reduced temperature 



L = Jk 

g\ 6L ' 



( 2 - 12 ) 



The 2+1 dimensional model has a deconfining second order phase transition at [24] 

t = t c = 0.614(3). (2-13) 



§3. Some results: 



Our main observable will be the correlation function of Polyakov loops L{x) and 
the corresponding screening mass i.e. inverse correlation length. 
In STALG we have 

L(x) = ^TrV{x) (3-1) 
V(x) = eW*)/^, (3-2) 

so that L(x) is a static operator and can be taken out of the integration over non- 
static variables. Therefore the comparison with the full model is straightforward. 

In Figure 3 we plot the screening mass as defined above. For details, see [17]. 
The screening masses in the naive tree level reduction do not agree with those of 
the full model. The agreement is, however, very good in the one-loop approximation 
form down to T ~ 1.5T C . At T c the screening mass goes to zero in the full model. 

The reduced model does not have this second order phase transition, and thus 
there will be necessarily a discrepancy in the values of the screening masses around 
T = T C . 

We have also measured the spacelike string tension. Also here there is good 
agreement between the full and the reduced model. In this case at high temperature 
there is also good agreement with the analytic result of pure gauge theory in 2d. 

In the 2d model one may in fact identify two screening masses, the lowest states 
m s of Tr<j) 2 (R = 1) and m ps of Trcj) 5 (R = —1). The corresponding correlation 
functions are nearly straight exponentials, corresponding to well isolated poles. The 
ratio m ps /m s between the two masses are between 1.5 and 2 depending on the 
temperature. Only a more detailed analysis could tell if there is agreement with the 
naive counting rule m ps /m s = 1.5 [25] For further details see [18]. 

The reduced model has two phases, a symmetric phase and a phase where the R- 
symmetry is spontaneously broken (Higgs phase). The transition in between is first 
order. The parameters in the one-loop approximation are in the unphysical broken 
phase, but very near the transition. In fact, we have performed the measurements 
at those values, but in the metastable symmetric phase region [17]. 
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Fig. 3. Physical screening masses Ms in units of g$ yT versus T c /T, in (2+1) D (black points) and 
2D (squares) . Also shown are the masses obtained with the tree level reduced action (diamonds) . 



§4. ^-symmetric dimensional reduction 

To obtain a reduced Z3 symmetric model we define instead of U, <j> effective 
variables U(x,i) and V(x) which are SU(3) matrices, where ^TrV(x) is the Polyakov 
loop. 

Define 



Z eff 

Su 
Su,v 
Sv 



S, 



eff 



d[U]d[V}e~ s ^f 

^ V ( 1 - -ReTrU(x, i)V(x + ai)U{ 
\ 2 Y,\TrV(x)/3\ 2 

X 

Su + Sjjy + Sv 



x,i) 



■V(x)- 



(4-1) 



The term Su,v is constructed such that developing V(x) in eft we get in lowest order 
the kinetic term of S2L above. Similarly, developing Sy we obtain terms proportional 
to Trcp 2 and higher orders. 

We investigated this model numerically, keeping A2 a free parameter [20, 21]. 
Similar models have been proposed in Refs [26, 27] and for imaginary chemical 
potential in Refs [28, 29]. 

At sufficiently high temperature in the deconfined phase this model should be 
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similar to the perturbative reduction. It has, however, an explicit ^-symmetry 
in the action, which may be broken at high temperature, but becomes restored at 
lower T. We have performed numerical simulations of the model, keeping Lq = 4, 
and scanning the (^3,^2) plane. With Lq fixed, the temperature is proportional to 
/?3. To characterize the phases we use the distribution of the Polyakov loop. 

I 1 1 1 1 1 1 1 n 1 1 




Fig. 4. A tentative phase diagram. Open points denote the location of peaks in susceptibility 
observed on 16 x 16 lattices. Black points denote the position of phase transitions established 
on 32 x 32 lattices with big statistics. The lines are hand drawn sketches of where transitions 
take place. 

We find a quite non-trivial phase structure, as shown in Figure 4. For A2 = 
there is in fact a transition to a confining phase at flz = 20, corresponding to r c = 
0.83, about 25% higher than in the full model. We checked that the screening mass 
actually vanishes there, as it should. Then fixing at 42 (r 2r c ), we made a 
more detailed investigation varying A2. The results for the susceptibility are plotted 
in Fig. 5. 

We find in fact two transitions, corresponding to three different phases. In 
the phase II of Fig. 4 with small A2, large the values of L are complex and 
concentrated near (1, 0) up to Z3 rotations, i.e. their phases are close to Himfi. This 
is the normal plasma phase. In the middle phase (phase III) L has also non vanishing 
values near the border of phase space, but with phases close to 2in(n + l/2)/3. This 
is a new phase. Finally for large A2, L « 0, as it is in the deconfined phase, and we 
propose that this phase (lb) is connected to the low temperature phase (la). 

It would be interesting to investigate, if the new phase plays a role in finite 
density QCD. 

We have investigated the behaviour of the screening mass. It is shown in Figure 
6. Since here A2 is a free parameter, we cannot predict the screening mass from the 
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x 30 




Fig. 5. Polyakov loop susceptibility as a function of A 2 for f3 3 — 42 on lattices of various sizes. The 
continuous lines were obtained through Swendsen-Ferrenberg reweightings. 



reduced model. If we demand on the other hand that the screening mass in the full 
and reduced models should be the same, we find this to happen at A2 ~ 0.18, well 
inside the physical phase II, in contrast to the earlier method. 
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Fig. 6. The screening masses ms for #3 = 42 and various lattice sizes. 



We have also supplemented the numerical investigation with an analytic mean 
field calculation, which gives qualitatively the same results, in particular the phase 
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structure. This calculation can easily be generalized to the 3+1 dimensional case. 
Again one observes the same phase structure. For further details on the Z3 symmetric 
reduction see [20, 21]. 

§5. Conclusions 

The usual method may be called perturbative reduction, although the reduced 
theory is calculated non perturbatively. In this method there is a systematic expan- 
sion in g 2 /T and where R is the typical distance of interest. No free parameters 
are introduced. As we have shown, for finite temperature SU(3) gauge theory in 2+1 
dimension this reduction works very well for screening masses and spacelike string 
tension down to T ~ 1.5T C . A difficulty of principle is that the parameters deter- 
mined from the reduction are not in the stable physical phase of the two dimensional 
model but in the metastable physical phase slightly beyond the transition. Fur- 
thermore in this method, as Z3 symmetry is explicitly broken, there is no confining 
transition. 

The Z3 symmetric reduction, which we propose, works all the way down to the 
transition, and the parameters can be chosen in the physical phase. 

However, in this case we have not fixed the parameters of the reduced model 
from the full model. In order to keep the quantitative predictive power, this has to 
be done. 

Qualitatively the reduced model is in good agreement with the full model. We 
have also found a Z3 symmetric new phase, which may be of importance at finite 
density. 
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